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Abstract 
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eral mechanical oscillators, coupled to the acoustic field in the fluid, is 
analyzed for different settings of the oscillators array. The dynamical 
models are formulated in terms of singular perturbations of the decou- 
pled dynamics of the acoustic field and the mechanical oscillators. 
Detailed spectral properties of the generators of the dynamics are given 
for each model we consider. In the case of a periodic array of mechan- 
ical oscillators it is shown that the energy spectrum presents a band 
structure. 
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Introduction 



In this paper we study the dynamics of a system consisting of one or more 
mechanical oscillators (the sources) coupled with the acoustic field they pro- 
duce in the compressible fluid surrounding them. 

Classical electromagnetism is perhaps the most well known case in theoret- 
ical physics where all attempts to construct a complete, covariant, causal, 
divergence free theory for the evolution of the fields together with their 
sources were unsuccessful up to now (in fact it is hard to say that there is 
a single case in Classical or in Quantum Physics in which this problem was 
completely solved). 

Whereas theories with extended rigid charges are quite well understood both 
at the classical and the quantum level (see e.g. the recent book for a 
systematic introduction to the subject and for a long list of references), there 
is no mathematically consistent theory of point charges interacting with their 
own electromagnetic field. Indeed Newton equations with Lorentz force 
require the fields to be evaluated at the particle positions, and this produces 
infinities due to the presence of the point-like sources. These difficulties 
directly lead to the need of mass renormalization. In his seminal paper Dirac 
[D] (also see |IWj . [K] . |M]). without using Lorentz force but exploiting the 
conservation of energy and momentum and considering their fiow through a 
thin tube of radius r, derived an equation for the motion of a charged point 
particle (the Lorentz-Dirac equation). As Dirac himself pointed out the 
equation obtained in the limit r | 0, together with the mass renormalization, 
leads to the presence of runaway solutions, i.e. solutions for which the 
acceleration increases beyond any bound even in the absence of external 
fields. 

An approach based on the theory of singular perturbations of the free dy- 
namics was initiated in |NPT| and [NP2] for the case of classical electro- 
dynamics of a point particle in the dipole (or linearized) case. Here the 
generator of the limit dynamics of both the field and the particle appears 
to be a singular perturbation of the generator of the free dynamics. The 
phenomenological mass plays the role of the parameter describing a suit- 
able family of self-adjoint extensions and the boundary condition naturally 
appearing in the domain of the generator results to be nothing else that a 
regularized (and linearized) version of the usual velocity-momentum relation 
in the presence of an electromagnetic field. In this framework runaway solu- 
tions are unavoidable because a negative eigenvalue appears in the spectrum 
of the generator after mass renormalization. 

Our interest in a similar problem in acoustics was prompted by the appear- 
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ance in 1999 of a paper by J. D. Templin |Tj. In that paper the author 
analyzed the dynamics of a simple model of a spherical oscillator interacting 
with the acoustic field it generates. The existence of a spherically symmetric 
radiation field (the acoustic monopole) makes the acoustic case significantly 
different from the electromagnetic one. Moreover the pressure field at the 
surface of the sphere completely characterizes the contact forces responsible 
of the interaction between source and field in the acoustic case. 
Templin performed a detailed analysis of the field emitted by the acous- 
tic monopole, explicitly computing both its radiation and near-field com- 
ponents. He noticed that a deduction of the reaction field obtained from 
the emitted radiation power, therefore neglecting the near field component, 
brings to an equation for the radius of the oscillating sphere showing run- 
away solutions. 

In analogy with what was done for the electromagnetic case in |NPlj we want 
to provide a formalization of the problem of a finite or infinite number of 
oscillators coupled with their acoustic field in terms of singular perturbations 
of the generator of the free dynamics. 

In this paper we will consider only the one dimensional case. In the abstract 
setting we will work in the physical model of interaction between sources 
and field will appear as the only possible extension of the free dynamics. 
The generalization to three dimensions is not straightforward. From one 
side a model of a physically relevant, symmetric, mechanical oscillator with 
finite degrees of freedom is lacking. On the other side point perturbations 
of the free dynamics are much more singular in higher dimensions. We plan 
to discuss the three dimensional case in further work. 

We want to stress an aspect of the dynamical system we analyze here which 
was extensively studied in different contexts. As an immediate consequence 
of the third Newton's law and of the assumption of persistent contact be- 
tween the fluid and the surface of the oscillators, the total energy, sum of 
the (positive) energy of the acoustic field Eac and the (positive) energy of 
the oscillators -Eoso is a constant of motion. As an immediate consequence 
one can exclude the existence of runaway solutions in this case. Moreover, 
lacking a mechanism of reflection of the acoustic waves at some exterior 
boundary, the motion of the oscillators should be damped and the energy 
should finally diffuse over the field degrees of freedom, for almost every ini- 
tial condition. The situation is reminiscent of the one investigated in |SWlj , 
[SW2. and |SW3j about the diffusion of energy from bound states to con- 
tinuous states triggered by time dependent perturbations in quantum and 
classical systems. In our system there is no external potential the interaction 
being given by internal forces. 
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This paper is organized as follows. 

In section^we introduce a list of notation and we briefly recall the equations 
for the acoustic field. Afterwards we exemplify the problem of the interaction 
between the field and a source in the completely solvable case of a single 
wall attracted toward the origin by a linear restoring force. 
In section [21 we analyze the case of a finite number of sources in the frame- 
work of the possible extensions of the free dynamics outside the points where 
the sources are placed. 

In section IHl we generalize the construction to the case of infinitely many 
sources and study the case of sources periodically placed on the real line. 
We give detailed results on the characteristic band structure of the spectrum 
of the generator of the dynamics. 

To the best of our knowledge this kind of systems of oscillators coupled with 
the acoustic field was never proposed and solved. A remark on the band 
structure of a similar model is in |GSj . 

1 The acoustic monopole in one dimension 

We give a detailed description of our model in the simplest case of one 
oscillator coupled with the acoustic field. 

Consider an infinite pipe filled with a non viscous, compressible fluid. We 
suppose that there is no friction between the fluid and the pipe and we choose 
a coordinate system with the j;-axis parallel to the axis of the pipe. The 
mechanical oscillator is made up of a very thin wall of mass M positioned 
in the pipe perpendicularly to the axis in x = 0. The thin wall is connected 
to a spring of elastic constant K. We analyze only one dimensional cases, 
hence the acoustic field is described by the pressure field p{x, t) and the 
velocity field The motion of the mechanical oscillator is described 

through the position and the velocity of the thin wall. 
The field p{x,t) represents deviations of the pressure in the point x at time 
t with respect to an equilibrium pressure Pq. In the linearized acoustics 
regime the continuity equation, the Newton's second law and the adiabatic 
equation of state read 



where p{x,t) is the deviation of the density in the point x at time t with 
respect to the equilibrium density pQ and a is the velocity of sound in the 




dp dv dv dp 




fluid. 
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Then we have for p{x,t) and v{x,t) the following coupled differential equa- 
tions 

dp 2 (^'^ 9v 1 dp 
dt ^ dx ' dt Podx ' 

We consider only small oscillations of the thin wall around its equilibrium 
position a; = 0, we indicate with y{t) the displacement of the wall from its 
equilibrium position at time t and we suppose that the wall remains always 
in contact with the fluid 

v{y{t),t) = ^ Vt>0. (1.3) 

Notice that we consider a wall of zero thickness. We make the approximation 
v{y{t),t) ~ v{0,t) and condition (|1.3|) becomes 

^(0,t) = ^ Vt>0. (1.4) 

The equation of motion for the position of the thin wall y{t) is 

My{t) = -Ky{t) - S {p{0+ ,t) - p{0' ,t)) (1.5) 

where S is the area of the transverse section of the pipe and we made the 

approximation p{y^{t),t) ~ p{0^,t). 

The total energy of the system is given by 



Etot — Eac + Eosc (1-6) 



with 



C r-oo c r-oo 

Eac = -^ p{xfdx + -^ V{xfdx (1.7) 

Ease = -y^ + Yy ' (1-8) 

Eac is the energy stored in the acoustic field while Eqsc is the energy of the 
mechanical oscillator. 

As the system is isolated the energy is constant. The motion of the wall 
produces acoustic waves thus transferring continuously energy from the os- 
cillator to the acoustic field. One then expects that y{t) decreases to zero 
when t oo. 

In spite of being a simple exercise, the exact computation of the solution of 
problem H1.2I) . (|1.4I) . (|1.5() and, in turn, of the damping rate of the oscillations 
rarely appears in textbooks. 
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In the following we give the solution of the Cauchy problem of coupled 
ordinary and partial differential equations with time dependent boundary 
conditions 

( = -aVo|^ Vt > Vx G M\{0} 

^ = --^ Vt>OVxGM\{0} 
ot po ox 

m = - ^ (p(0+, t) - ^(0", t)) Vt > 
p[x, 0) = f{x) Vx E M\{0} (1-9) 

v{x,0)=g{x) VxGM\{0} 
2/(0) = yo 

y(o) = m 

v{0,t)=y{t) Vt>0 

where f{x) and g(x) are two real functions and loq = K/M. 
Suppose that 



f{x) e C^{R) ; g{x) G ^^(M) and yo = ^^/'(O) ; Vo = 9{0) , (1.10) 



1 

then the solution of problem (|1.9() reads 



p{x, t) = pf{x, t) + apo sgn(x)y It ) (1.11 



a 



v{x,t) = Vf{x,t) + Y [t-^-^) (1.12) 



yit)^-^-I^J-^ (1.13) 



where pf{x,t) and Vf{x,t) are solution of the wave equation in (— C)0,+oo) 
with initial conditions f{x) and g{x) 

(1.14) 

, , gix — at) + gix + at) 1 , . ^, 

Vf{x, t) = ^ ^ + ^ (/(x - at) - fix + at)) , 

F{t) = -VfiO,t)-u;^oVfi0^t) (1-15) 
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and 

f F{t')e''*'-'-'"'dt' - f Fit')e'''"-''1dt' 
= 

with (5± = (-7 ± A/72 -4a;^)/2 with 7 = 2apQS/M. 

By conditions (|1.1U|) one easily obtain that y{t) and y{t) are both continuous 
and decrease exponentially to zero with decay constant r = 7/2. 

2 Singular perturbations of the free dynamics 

In this section we present a generalization of problem (|1.9jl formulated in 
terms of a unitary flow on a space of finite energy. 

Let us consider a system of n thin walls positioned in a pipe, perpendicularly 
to its axis. Let S = {si, . . . , s^} C M be the set of equilibrium positions 
of the thin walls. The i-ih thin wall, placed in Sj, has mass Mj and is 
connected to a spring of elastic constant Ki . The acoustic field is described 
by the pressure field p and the velocity field v. The motion of the walls is 
described by the displacements yj from their equilibrium positions and by 
the corresponding velocities Zj. 

The generator of the dynamics, A, will be defined as a singular perturbation 
of the skew-adjoint operator A generating the uncoupled evolution of the 
acoustic field and of the oscillators. 
The system of first order differential equations 

| = -<.V.| V«M (2.1) 



dyj 

dt ^ 
dzj_ ^ _Kj_ 

dt Mj 



1 < j < n (2.3) 
yj l<j<n (2.4) 



describes, in the linear approximation, the independent evolution of n me- 
chanical oscillators and of the acoustic field. 

We want to show first how equations (|2.1|) - (|2.4|) define an unitary fiow in 
a complex Hilbert space. To this aim let us consider equations (|2.1() - (|2.4() 
for complex functions v, p, yj, zj of position and time. 
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The set of all the displacements and velocities will be represented respec- 
tively by the vectors in C" 

y_ = yiCi H h yne„ , z = zie^ H h 2;„e„ , (2.5) 

where 6;^, . . . , e„ is the canonical orthonormal basis in C". 

Let us denote by (M) the space of square- integrable functions on the real 

line. H^(W) indicates the homogeneous Sobolev space of locally square- 

integrable functions with square-integrable (distributional) derivative, and 

ffi(M) the usual Sobolev space H'^(R) := H'^{R) n L'^{R). 

Therefore the hnear operator A in L^(R) L'^{R) ® C" © generating the 

dynamics (|2.1|) - 1)2. 4|) is defined by 

A : H'^{R) © H^{R) © C" © C" ^ L^{R) © L^{R) © C" © C" (2.6) 




A{p,v,y,z):=\-a'po^,-^^^,z,- }] ^ y, e, | , (2.7) 

where a, po, Kj, Mj, 1 < j < n, are the positive real constants representing 
the physical parameters. 

In the following a capital Greek letter will indicate a generic vector {p,v,y,z) 
inL2(K)©L2(M)©C"©C". 

^ is a real operator, i.e. it preserves the (physical) linear subspace of real 
elements 

{iP,v,y,z) ■■ p{x) G M, v{x) G M, y G M", z G M"} . (2.8) 
A is skew-symmetric with respect to the scalar product 

{{^1,^2)) = ^ {pi,P2)+Po{vi,V2) + - V Kjyijy2j+MjZijZ2j , (2.9) 
a Po -J -,7-^ 

l<j<n 

where (• , •) indicates the standard scalar product in L^(M), S is the area 
of the transverse section of the pipe and ~ denotes complex conjugation. 
L^{R) © L^{R) © C" © is a Hilbert space with the scalar product (OH . 
The square norm of a vector = ((^', ^')), defines the total energy of 

the system in the state \E' 

^^tot = |||^f =Kc + ^o.c (2.10) 
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where Eac is the energy stored in the acoustic field while Eosc is the energy 
of the oscillators 



l<i<n 



For any (" G C\'iM the resolvent of A is 

{-A + crHp,v,y,z) =\po( -^ + ^] ( -^ + ^P 



1 / d 



(2.11) 



dx"^ J \ dx a^pQ 
,eyW dp \ ^ M^z,+CM,y 



■2 



«Vo \ dx^ ^ , , „^ , ^^^^^^ 



l<j<n J ' ^ 1 



Since 

Ran(-^ ± 1) = © L'^{R) © C" © , (2.13) 

A is skew-adjoint. Moreover the essential spectrum of A is purely absolutely 
continuous and 



aessiA) = aaciA) = iR , app{A) = |±i , 1 < j < n| . (2.14) 

Being skew-adjoint the operator A describes, by Stone theorem, the un- 
coupled evolution of the acoustic field and of the oscillators through the 
unitary fiow exp tA corresponding to the Cauchy problem for the first order 
differential equation 

|^(t) = A^{t) (2.15) 

which is equivalent to the system written at the beginning of the section. 
Now we consider the linear operator obtained by restricting A on the set 
of vectors in its domain satisfying: 



{v{sj) = Zj, l<j<n} (2.16) 

which represents the kinematic constraint (|1.4|) at each thin wall. Aq is a 
closed, densely defined, skew-symmetric linear operator with defect indices 
(n,n). We want to characterize the skew-adjoint extensions of Aq. The 
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family of extensions of Aq can be parameterized by relations K C © 
which are skew-symmetric, i.e. such that K = (j^i^)-*-, where ^{z_i,z^2) '■— 
(l2)li) (see e.g. |GGj . Theorem 1.6, chapter 3, for the analogous self-adjoint 
case). A skew-symmetric relation in C" ® C" extends the notion of the 
graph of a skew-symmetric operator G : C" — > C" through the relation K = 
{{z,Qz) , z G C"}. In order to be a candidate to describe the interacting 
dynamics of the system under analysis a skew- adjoint extension of Aq must 
be local and real i.e. it must generate a coupling between the fields evaluated 
in Sj and the j-th oscillator, 1 < j < n, and it must preserve the linear space 
of physical data defined in (|2.8|) . The only admissible extension different 
from A itself will be the one corresponding to the graph of the zero operator, 
= 0. The next theorem completely characterizes such an extension: 

Theorem 2.1. The only local, real and skew- adjoint extension ofA^ is given 
by 

A : D{A) c L^{R) © L^{R) © C" © ^ L^{R) © L^{R) © C" © C" , (2.17) 



D{A) = = {p,v,y,z) : p £ L"^ {R) n H'^ {R\S) , v £ H'^ (R) , y G C", 
z G C^, p{sf) -pisi) = ai, v{sj) = Zj, a G C"} , 

(2.18) 

Mp^v,y,z) := 

I 2 dv I dpo ^ [Kj S \ \ (2.19) 

- [-^ « d^- -« d^- - (m- w + 5^ • 

Here po G H'^{R), 

Po{x) ■■= p{x) - ^ J2 o-j sgn(x - Sj) , (2.20) 

l<i<n 

denotes the regular part of p. The resolvent of A is given by 

{-A + cr' = {-A + c)''+ (r(c)-i),^. G^©GJ, CgCVm, 

l<i j<n 

(2.21) 

where 
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and 



(2.23) 



GUx) 



-G'Ax-Sj). 



Gl^ix - Sj) , 



Kj + C^Mj ' ivTj + CMj J ' 

(2.24) 



g^{x) = ±—e^^\''\'\ g^(x) = --sgn(x)e^^l"l/^ ±ReC>0. (2.25) 

Proof. Since iAq is a closed, densely defined, symmetric operator with defect 
indices (n, n), all its self-adjoint extensions can be obtained by the famed von 
Neumann theory on self-adjoint extensions (see e.g. Theorem X.2 in |RS2j l. 
However, since Aq is obtained by restricting the skew-adjoint operator A to 
the kernel of the continuous, surjective linear operator 

r : (M) (M) C" © C" ^ C" , (2.26) 



r{p,v,y,z) := ^(1^(5^)-^;^)%, (2.27) 

l<i<n 

it is easier to make use of the (equivalent) procedure developed in |Plj (also 
see the appendix in |P2j for a compact review). Here below we provide the 
(almost) self-contained construction of the skew-adjoint extensions of Aq by 
using such a procedure. 

Given the map r we can define the bounded linear operator 

G{C) := t{-A + C)~^ : L^R) © L^(R) © C" © C" ^ C" . (2.28) 
By the relation 

G{C){p,v,y,z)= {{^,{p,v,y,z)))ej, l<j<n, (2.29) 

l<i<n 

G{C) is represented by the vector G;^. By G'(C) we define the bounded linear 
operator 

GiO := - G(-C)* : ^ L\R) © L^{R) © C" © C" , (2.30) 
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where G{C)* indicates the adjoint of G{C)- The action of G{() on C" is given 
by 

GiOe^=Gl, l<j<n. (2.31) 

Let us notice that 

Ran(G(C)) n H\M.) © H\R) © C" © C" = {0} . (2.32) 

Now we consider the hnear operator r0(C) : C" C" represented by the 
matrix Qij + T((^)ij, where : C" — > is skew-Hermitian. By noticing 
that 

r0(C)-re(O = T(G(e)-G(C)) (2.33) 
and that, by the definition of G{C) and by the first resolvent identity, 

(C - C)i-A + Cr'GiC) = G{C) - G(C) , (2.34) 

one has that re(C) satisfies the identity 

r0(C) - re(0 = (C - OGiOGiC) ■ (2.35) 

By the definitions of G{C) and G(C), by ^IT^ and by re(C)* = -re(-C), 
it follows that detre(C) / for any C G C\iR and that 

RiC) := i-A + C)-^ + GiOTeiCr'GiC) (2.36) 

satisfies the first resolvent identity 

{c-omRic) = m-m (2.37) 

and 

R{0* = -R{-0 (2.38) 
(for details see |Plj ). Moreover R{C) is injective by (|2.32|) . Therefore 

A:= -R{C)-^ + C (2.39) 

is well defined on 

D{A) := Ran(i?(C)) . (2.40) 

By (|2.37)) such a definition of A is (^-independent. A is skew-symmetric by 
1)2. 38() and is skew-adjoint since 

Ran(-i ± 1) = l2(M) © L^{R) © C" © C" (2.41) 
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by construction. 

Since we require A to be real, i.e. to preserve the linear space ()2.8() . we 
have to restrict the choice of Q to real, skew-symmetric matrices. Off- 
diagonal elements in the matrix G would correspond to non local couplings 
between the pressure field and the oscillators. Since we are looking for local 
interactions the only admissible choice for the skew-symmetric matrix Q is 
9 = 0. 

By (ITini) (p, v,y,z) G D{A) if and only if 

p{x)=p^{x)- ^ (TiCr%{v^{sj)-z^,)g'^{x-s,), (2.42) 

l<i,j<n 

= v^{x) + ^ {nO-%{v^{s,)-z^,)gc{x-Si), (2.43) 

,=2,-s E (r(cr')«^g^... P.44) 

l<i,j<n 

l<i,j<n « S I 

with {p({x) ,ff(x) ,y^,z^) G D{A). Posing 

Mp^v,y,z) = {Ai{p,v,y,z),A2{p,v,y,z),A3{p,v,y,z),A4{p,v,y,z)) , 

(2.46) 

The action of A on {p,v,y,z) is given by 

[Ai{p,v,y,z)]{x) = - a^po^ix) - C ^ {T{Cy^)ij {v(;{sj) - z^j)g'^{x - s, 

l<i,j<n 

(2.47) 

[A,ip,v,y,z)]{^)=-^^i^) + 4- E im-'h (vdsj) - za)Od^ - 

(2.48) 

A3{p,v,y,z)=z^-CS (r(C)-^)^ /5f'f ~ y (2-49) 

l<i,j<n * ^ * 



(2.50) 
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By the definitions of D{A) and T(Q one has 

Ai {p, V, y, z) = -a^po ^ (2.51) 

Mp,v,y,z) = z (2.52) 

and, defining 

a, := p{st) - p{sT) = iX{C)-% {vd^,) - za) , (2.53) 
i<i<" 

formula H2.5U() becomes 

l<j<n J l<i<n ^ ' ^ * 

Then, posing 

p(x)=p^(j;)- ^ cFjG'd^-Sj) =po{x) + - ^ cTj sgn(x - Sj) , (2.56) 

l<j<n l<i<n 

one obtains 

[^2(.P,v,y,z)]{x) = 



(2.57) 
(2.58) 



1 f^Po 



Finally 

t;(^.)=t;c(^.) + ^ 5] (r(c)^^).,(^c(^.)-^c.)^e^^i^'=-^^i/'^ = 

(2.59) 

=vdsk)- E (r(C)-%-(^f(«.)-^c.) ((r(C)).. + ^^^^) = 

(2.60) 
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□ 



By the previous theorem the differential equation 

d 



^(t) = A^fit) (2.62) 



is equivalent to the system of equations 



9p 2 dv 



dv _ 1 dpo 
dt po dx 

dg _ 

Ht ~ - 




(2.64) 



(2.65) 



and the corresponding Cauchy problem generates the strongly continuous 
unitary group of evolution exptA on L'^(R)©L'^(M)©C"'©C" which preserves 
D[A). Here 5s ^ denotes the Dirac mass at the point Sj and aj (see ()2.18() 1 
is the pressure jump at Sj. 

It is worth noting that the only real, skew-adjoint extension of the free oper- 
ator A restricted to the space of the vectors (p, v, y,z) such that v{si,t) = Zi 
corresponds to the relevant physical coupling between the pressure field and 
the oscillators. 

The next result will be useful in the spectral analysis of A. 
Lemma 2.2. The matrix 

r±{X)-^ := limr(A ± e)-^ (2.67) 

is well defined for any A G iM\{0}. 

Proof. We give the proof only for the matrix r+(A). The proof for r_(A) is 
analogous. 

Let the matrix r_|_(^) be the analytic continuation to C\U"^]^ {zizi^yKj/Mj} 
of T{() defined for ReC > in (|2.22|) . Suppose that Si > Sj if z > j, then 

r+(C) = -n(c)-r(c) (2.68) 
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where 11 is the operator 

n= (r(C)®^+(C)) (2.69) 

with ^^(C) = Yli ^^2apo -i' ^^^^^ ^(C) is the upper triangular matrix 



ann — 



T{Qij = { Ki+c^Mi ^ apo " (2.70) 
[0 i> j 

We use the formula 

' = ~n(C) + T(c) ^~W)^W) ^^^^ n(c) + T(c) = ^^''^^^ 

-1)" / 1 
^(^n(C)^ 



n=0 



valid for all C, for which the series converges. 

Matrix T{C,) is invertible and its inverse T{C,)~^ is a lower triangular matrix 
with (r(C)-^),^ = 1/ (T(C)),,. The eigenvalues of r(C)-^ are 1/ (r(C)),, and 
we can write 

T(c)-^ = D(c)r(C)-^D(C)-' (2.73) 

where -D(C) is a unitary matrix, analytic for C, G C\{0} and 
We obtain for r+((^)'~^ the expression 

oo 

r+(C)-^ = -D{0 E(-l)" {±{0®x{QynO-'D{Q-^ (2.75) 

n=0 

with 

too), = ^i±|^ (i^(C)-V-(O), (2.76) 

(X(C)), = (I)(C)- V(C)), (2.77) 

Then 

^ oo 

r+(C)-^ = -r7n+E(-i)" ((x(C),i(C))c")"i)(C)V:(C)c5x(C)T(C)-'i?(C)-^ 

n=0 

(2.78) 
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For all C for which the series converges one has 



Consider the scalar product in C" 

ixiOMO)^. = E (-P(C)-^^+(C)), ' ' {DiO-'riO), ■ (2.80) 

i=l ^ 
Notice that, for A G iM\{0}, (x(A), V'(A))^„ G and 



-i(x(A),V^(A)>^„ ^+cx) forA^+icx) (2.81) 
-i(x(A),V:(A))^„ ^-oo forA^iO+. (2.82) 

Then there exists at least one point A G iM in which (x('^)) V'(-^))c" = 0. 
In a neighborhood of this point the series converges and defines an analytic 
function. By (|T7^ and (jTHni) it is clear that r+(C)"^ exists for any C G 
C\{0}. The same relations show that one can put r+(C)~^ := if C = 
i^Kj/Mj,j = l,...,n. □ 

The following theorem completely characterizes the spectrum of A. 

Theorem 2.3. The essential spectrum of A is purely absolutely continuous 
and 

aess{A) = aac{A) = iR , app{A) = {0} . (2.83) 
Any vector of the kind 

\l ^Sn(x - Sj), 0, - E ^J^i' ^ ' (2-^^) 

\ i<i<n i<i<" ■' J 

with 

is an eigenvector corresponding to the (n — l)-fold degenerate eigenvalue 
A = 0. 

The generalized eigenf unctions i^(A) corresponding to the point of the ab- 
solutely continuous spectrum relative to right (+) and left (— ) incidence are 
given by 

$±(A,x) = (0±(A,x),0±(A,x),0j(A),0f(A)) AGiM (2.86) 
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(2.87) 



l<jj'<n 

A ST^ p±\sj/a 
l<i,j<n 



with C = y/apo/{4:TT). 

Proof. For C G p(^) n {-A + C)"^ - (-^ + C)~^ is of finite rank, then 

from Weyl's criterion (see e.g. |RS4j Theorem XIII. 14) one has cJess(^) = 
(Tess(^) = iM. Moreover, by Birman-Kato invariance principle, the wave 
operators ^1±{A, A) exist and are complete (see e.g. |RS3j . Corollary 2 to 
Theorem XI. 11). Thus aac{A) = aac{A). 

Let fi^ be the singular continuous part of the spectral measure on iM. cor- 
responding to A and Since ||G(C)^|| < oo for all C £ C\app{A) and for 
aU * G 

D := {"^ = {p,v,y,z) : p e L'^ (R) H L"^ (R) , v e (R) H (R)} , 

(2.91) 

by Lemma lO and |RS4j . Theorem XIII. 19, one has supp/if' C {0}Llapp{A) 
i.e. supp fi^ = since fi^ has no atoms by its definition. Since D is dense 
this gives (7sc{A) = 0. 

One can check that any vector ^ of the kind (|2.84|) is in the domain of A 
and solves the equation A^ = 0. The degeneration of eigenvalue {0} follows 
from condition l|2.85jl . 

Suppose now A G iM\{0} and consider the equation A^ = A^. This pro- 
duces, if = {p,v,y,z}, the equation 



v"-^v 



A2 A 



i<i<n 



E ^i-^^. ' (2.92) 



with fjj G C, « = 1, . . . , n, which has no square integrable solution. 

The expression for the generalized eigenfunctions is a consequence of the 

Stone's formula (see e.g. |RSlj . Theorem VII. 13) which gives the generalized 
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expansion formula 



1 . 

* = s- lim s-lim— / [R{X + e) - R{\ - e)]"^ dX . (2.93) 

aj,-oo,fetoo ej,0 iVT 

□ 

In the following lemma the asymptotic behavior of the oscillations of the thin 
walls is characterized. It is proved that the oscillators relax (as \t\ oo) 
towards their equilibrium positions for any initial data orthogonal to the 
eigenspace relative to eigenvalue zero. For example this is true for any 
initial datum of the kind ^'o = ip,v,0,z) where the support of p is outside 
the interval containing the points si,. . . ,Sn which denote the equilibrium 
position of the walls. 

Lemma 2.4. Given orthogonal to the eigenspace relative to eigenvalue 
zero, let us denote by {y{t),z{t)) the projection onto C^SC" ofe^'^^Q. Then 

lim ||2/(t)||c" = c-i^d lim ||l(i)||c" = . 

Proof. Let P{dk) be the projection-valued measure corresponding to the 
self-adjoint operator —iA. Since \I'o is in the absolutely continuous sub- 
space, for any ^' the bounded complex measure ((^', P{dk)^o)) is absolutely 
continuous with respect to Lebesgue measure and hence its density belongs 
to L^(M). Thus, by the spectral theorem and Riemann-Lebesgue lemma, 

lim ((^',e*^*o)) = lim / e''^^ {{^ , P{dk)^o)) = . (2.94) 

\t\-^oo |i|^oo Jr 

By taking ^' = (0, 0,6^,0) and = (0, 0,0,6^), i = 1, . . . , n, one then obtains 
lim yi{t) = and lim Zi{t) = . (2.95) 

|t|-^oo |t|— +00 

□ 

In order to obtain more precise estimate on the asymptotic behavior of solu- 



tions of equation (|2.62|) . for particular initial conditions, a detailed analysis 
of r(A)~^ is required. For example in specific cases one can prove existence 
of frequencies which are totally transmitted by the array of oscillators. 
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3 Kronig-Penney model in acoustics 

It is possible to extend the previous construction to the case of an array of 
infinitely many oscillators. We prove that in the case of a periodic array of 
identical oscillators the energy spectrum shows a band structure. 
As a first step we define the operator A introduced in section |21 when S = 
{si, S2, . . . } is a denumerable set such that 

d := inf |si - Sjl > j G N. (3.1) 
Defining the linear map 

oo 

T{p,v,y,z) := ^(w(sj) - Zj)ej , (3.2) 

where {gj}^ is the usual complete orthonormal system for P, one has the 
following 

Lemma 3.1. r is bounded as a map on H'^{R) © F^(IR) ®P ®P to P. 

Proof. We will follow closely |A(;H-KHj . Let {Ij}"^ be a partition of M and 
let K(x — y) be the kernel of (—A + 1)~^/^. Since 

oo „ 

v{x) = V / K{x - y) [(-A + lfl''v]{y) dy , (3.3) 

to prove the lemma amounts to show that the infinite matrix 

\ 1/2 

Mij := [ I Kix- Sif dx (3.4) 



corresponds to a bounded linear operator M on By Lemma C.3 in 
|A(;H-KH| . one has 

oo / \ 1/2 oo / \l/2 

||M||p_;2 <sup^ f y K{x-Sifdx\ swp'^ij K{x-Sifdx\ 

(3.5) 

Since 

du 

^ x>0, 



y/x vr 7o x + fji 
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by functional calculus one has 

1 /"oo p-\J^+p? \x-y\ 
K{x-y) = -\ ^ dii. (3.6) 

By taking Ij = [sj — ej , sj + 5.j ) , where ej is one half the distance between 
Sj and the preceding point and 6j is one half the distance between sj and 
the successive point, a straightforward calculation leads to 

/ K{x-Sifdx<-^e-\'^-'^\/^, (3.7) 
Jij d^T^^ 

from where the estimate ||M||p p < +oo follows immediately. □ 

The construction proceeds now along the same lines as in the case of a finite 
set of points. We state the final result: 

Theorem 3.2. Let {Kj}"^ {Mj}^ , Kj > 0, Mj > be in 1°^ and suppose 
that {Kj/Mj}f and {\/Mj}^ are in l°° too. The linear operator 

A : D{A) C L^{R) © l2(R) ^2 ^2 ^ ^2^^^ ^ ^2^^^ Ql'^Ql'^, (3.8) 



D{A) = {{p,v,y,z) : peL^{R)nH^{R\S),veH^{R), yel^,zel'^, 
P{4) - Pi^i) = '"(sj) =Zj, ae , 



(3.9) 



^ dv 1 dpo ^fKj S \ \ (3.10) 

is real and skew- adjoint. Here po G H^{R), 

1 °° 

Po{x) := p{x) - 2 XI ^Sii*^^ ~ ' ^^■^^'> 
denotes the regular part of p. The resolvent of A is given by 

oo 

(-i+C)-^ = (-^+C)-'+ X (r(C)-'),,- Gi^Gi, c e CVM. (3.12) 
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Now we can proceed to the study of a periodic system. We use the same 

notation of gGH-KH . 

In this case S will be the "Bravais" lattice, 

S = {nL:n£Z}, L>0, (3.13) 

and S the "Brillouin" zone, 

We consider a Hilbert space W on © © © P in which the scalar 
product is defined by 

1 KM 

-^{pi,P2) + po{vi,V2) + -w{y.,yn) + -7r(li>l2> (3-15) 

where (•, •) represents either the usual scalar product in L^, when concerning 
pressure and velocity fields, or the usual scalar product in P, for y and z. 
M, K and S are positive constants representing the mass of oscillating walls, 
the elastic constant of the springs and the area of the transverse section of 
the pipe. 

The Hilbert space Ti can be decomposed as 
n = W-^n{S, b-^dO; l2([-L/2, L/2)) © l2([-L/2, L/2)) © C © C) 

(3.16) 

r® dft f \ 

= W-^ — L\[-L/2, L/2)) © l2([-L/2, L/2)) © C © C 

J[-b/2,b/2) b y J 

(3.17) 

where 

#: n ^n{S,b-'^ de;L\[- L/2, L/2))® L'^{[-L/2, L/2))® C®C) (3.18) 
^{p,v,y,z)= ({Wp){d,u),{Wv){e,u),{Wy){d),{Wz){d)) (3.19) 
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iWp){0, v) = p{e, i^) = Y, e"%(z^ + nL) (3.20) 

{Wv){9, u) = v{9, iy) = Y^ e*"^^w(i/ + nL) (3.21) 

(%)(^)^y(0) = ^e-%„ (3.22) 

iWz){0) = z{9) = e^"^^z„ 1/ e [-L/2, L/2) , ^ G [-6/2, 6/2) 

(3.23) 

and 

: n{S,b-^de; L'^{[-L/2, L/2))(BL'^i[-L/2, L/2))eC®C) ^ n (3.24) 



^ ^(p,^',y,I) = 

= ((W^-ip)(z. + nL), {W-^v){iy + nL), {(W^-iy)„}, {{W-^z)n} 

(3.25) 

{W-^p){u + nL) = b-^ dee-'''^^p{e,v) (3.26) 
{W'^v){u + nL) = b-^ dee-''^^^v{e,u) (3.27) 



(VF-iy)„ = b-^ / dee-'''^^y{e) (3.28) 

7-6/2 
/.fe/2 

{W-^~z)n = b-^ dOe-'^'^^zie) e[-L/2,L/2) , neZ. (3.29) 

The scalar product in L'^{[-L/2, L/2)) ® L'^{[-L/2, L/2))®C®C is defined 
by 

1 , , , ^ K - M - , , 

(P1,P2)l/2 +P0(?^1,';^2)l/2 + -Wmm + -^ZiZ2 (3.30) 

where (•,-)l/2 indicates the usual scalar product in L^([— L/2, L/2)). 
From Theorem 13.21 we obtain the following 

Corollary 3.3. The linear operator 

A: D{A) cn (3.31) 
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p{nL'^) — p{nL ) = cr„ , v{nL) = Zn^n , £ € 



(3.32) 



2 dv 1 dpo K S 



A{p,v,y,z) := [-a Po^,--^,z,--y- -a) , (3.33) 
where the regular part ofp{x), denoted with po G 5'^(M), is 

Po{x) =p{x) an sgn{x -nL) , (3.34) 
nez 

is real and skew- adjoint. 

We want to study the spectral structure of A. To this aim we introduce the 
family of operators A{9) 

A{e) : D{A{e)) C L2((-L/2, L/2)) © L^{{-L/2, L/2)) © C © C 

o o (3.35) 

^ L^{{-L/2,L/2))®L^{{-L/2,L/2))®C®C 



D{A{e)) = i^{p{e),m,m,m) ■ 

m G H\{-L/2,L/2)\m, m e H\{-L/2,L/2)), 
y{e) G C, ~z{e) G C 

p{e,^+)~p{9,Q-)=a{e), 
v{e,o) = z{e), a{e)ec, 



b b 
2' 2 



(3.36) 



2 d5(9) 1 dpo{9) 5 \ (3.37) 
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where po{6) G H^(M.) is the regular part of p{6) 

Po{e, u) = p{e, v) - h{e) sgn(i.) . (3.38) 

Boundary conditions for p{6, v) and v{6, z^) in = and u = ibL/2 are such 
that ah operators in this family are skew-adjoint with respect to the scalar 
product (121301) • 

The operator A is related to A{6) by the relation (see |A(tH-KH] ) 

WAW-^ = / l-^(^) • (3-39) 

J[-b/2,b/2) 

The spectrum of A{6) is described by the following 

Theorem 3.4. Let 9 £ [—6/2,6/2) then the spectrum of A(9) is purely 
discrete, in particular its eigenvalues En{9) are given by 

Er,{9) = K{9) = 2iU0)j ; n E Z , en(e) e M (3.40) 

where £,n{9) are the real solutions of 

9L 9L 
sin ^ [sin ^ — F(C) cos ^] cos^ — = cos ^[cos^ + F(^)sin^]sin2 — (3.41) 

The corresponding eigenf unctions are 

^n{e,x) = {pn{e,v),Vn{9,v),yn{e),Zn{0)); n G Z, € [-6/2, 6/2) 

(3.43) 
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PniO^f) =Cn 



sin Un - ^ - F{e,n) COS Kn - ^ COS 



2 J 



9L\ 



2 J 



L 



i sin ( ^1 

iCn 
apo 



2Cn 



2) [-^^—--nCn)-COS—U 



(3.44) 



2 / 



sin ( - ^ ) - F{^_ri) COS ( C„ - ^ ) ) sin '^^+ 



2 J 



+ i sin ^,^„ 

. CnL 



9L\ 



2 y 



2e 



cos ——1/ + sin 



L 



L ' 



Zn{0) = sin Un - — 

apo V 2 y 



sin (Cn - ^) 
0L\ 



(3.45) 

(3.46) 
(3.47) 

(3.48) 



For 9 G [— &/2, &/2) zero is an eigenvalue with eigenfunction 

[ ^ f OL 9L ^ ^\ ^ S 9L 

*o = ( Co (cos — - ism— sgn(i^) I , 0, 2?Co— sm — , 

Moreover the following chain of inequalities holds 

< Ei{0) < Ei{-b/2) < E2{-b/2) < ^2(0) < ^3(0) < E^[-h/2) < 
< Ei{-h/2) < ^4(0) < -£5(0) < Ei{-h/2) < E^{-b/2) <... 

(3.49) 

In general the eigenvalues En{9) are all distinct and non degenerate. If 
^o/^g = ^/2 with n G N there is just one two fold degenerate eigenvalue 
equal to mr/2, such eigenvalue corresponds to 9 = for n even and to 
\9\ = 6/2 for n odd. 

If E{9) is an eigenvalue then —E{9) is an eigenvalue. 
Given 9 G [—6/2,6/2) the following relation holds 



En{-9) = En{9) . 



(3.50) 



Proof. Eigenvalues and eigenfunctions (|3.4()|) - (|3.48|) are given by direct com- 
putation. We solve the system of equations 



Ai9m9),v{9),m,m) = x{p{9),v{9),m,m) a g 



(3.51) 
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with the condition v{6,0) = z{9), the solution reads 



P{e,u) =C(e) cos ^ + 2^(0 

via, v) =- sm — : 

tapQ L lapo 



sm — i^(^)sgn(z/j cos — — 

Li Li 



COS — — h r [t) sm — - — 

±J ±J 



(3.52) 
(3.53) 



where ^ = —iLX/{2a) G M, C(.^) and -D(^) are two unknown functions of 
^. To determine C(^) and D(^) we have to take into account the boundary 
conditions 



(3.54) 



This system has only the trivial solution C{(,) =0 and D{S^) = for the 
values of ^ for which the determinant of the matrix of the coefficients of the 
system is zero. The condition that the determinant is zero implies equation 
1)3.41(1 for the eigenvalues. For ^ satisfying condition (|3.41() the solutions of 
the system of dependent equations 1)3. 54() give the eigenfunctions. 
For 9 = and 6 = —6/2 relation (|3.41|) becomes 



tan^ = or 
cot C = or 



tan^ = F(^) 



cote = F(0 



M 
M 



0^1 
21 



IT 



(3.55) 
-6/2 (3.56) 



Graphic solutions of the transcendental equations ()3.55|) and ()3.56|1 are given 

The chain of inequalities ()3.49() 



in the upper part of figures 1(a) 



and 



1(b) 



follows by the monotone behavior of F{^). 

Degeneration of eigenvalues for coo/^g = the fact that —E{6) is an 
eigenvalue if E[9) is an eigenvalue and relation ()3.5fl|) follow directly by 
equation (jTm and by F(^) = -F(-^). □ 



One can show that there is a band structure writing equation ()3.41|) as 



2 0L 
tan — = tan^ 



tan^-F(0 



1 + F(^)tan^ 



(3.57) 



It is possible to find solutions of equation H3.57() only for values of ^ such 



that the r.h.s. is positive. In the lower part of figures 1(b) and 1(a) the 
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(a) With non degenerate eigenvalues. 

M/Mg = 0.5, UJo/Ug = 1.2 




Figure 1: The upper part of figures shows the graphical solution of equations 
H3.55() and H3.56() . The lower part shows the band structure. 



resulting band structure is shown. The figures clearly show that the width 
of the gaps is connected to the structure of the spectrum. In particular 
figure 1(b) shows that when there is a degenerate eigenvalue, oJo/ojg = n-n /2 
with n S N, a gap disappears because of the overlapping of two bands. 
The bandwidth increases, when the ratio M/Mg decreases. 
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